IPPP/10/83; DCPT/10/166 



Lower Bounds on the Smallest Lepton Mixing Angle 

Shamayita Ray, ^'Q Werner Rodejohann,^']^ and Michael A. Schmidt'^'[^ 

^Laboratory for Elementary-Particle Physics (LEPP), Cornell University, Ithaca, NY 14853 USA 
^ Max-Planck-Institut fur Kernphysik, Postfach 103980, D-69029 Heidelberg, Germany 
"'institute for Particle Physics Phenomenology (IPPP), Durham University, Durham DHl 3LE, UK 

We give minimal values for the smallest lepton mixing parameter |?7e3j, applying 2-loop renormal- 
ization group equations in an effective theory approach. This is relevant in scenarios that predict an 
inverted neutrino mass spectrum with the smallest mass and |?7e3| being zero at tree level, a situa- 
tion known to be preserved at 1-loop order. At 2-loop, |f/e3| is generated at a level of 10~^^-10~^*. 
Such small values are of interest in supernova physics. Corresponding limits for the normal mass 
ordering are several orders of magnitude larger. Our results show that | Ue3 \ can in general expected 
be to be non-zero. 
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The subject of lepton mixing and neutrino mass has 
entered the precision era. It is confirmed that there are 
three flavors of neutrinos which mix among themselves 
because the flavor states are not identical to the states 
with distinct masses mi, 7712 and 7713. The parameters de- 
scribing the mixing are the three angles 9i2, ^23, ^13 and 
a phase S. Two more phases (^1,2 are needed if neutri- 
nos are Majorana particles, which is the case in basically 
all extensions of the Standard Model, and also from an 
effective field theory point of view. The two mixing an- 
gles 012 and ^23 as well as the mass squared difference 



- mf and the magnitude of Am 



have been determined with high accuracy, and further 
improvement is foreseen with a variety of experimental 
approaches l\. However, of particular interest is the 
parameter \Ue3\ ~ sin^ia, which describes the electron 
neutrino content in the heaviest (lightest) neutrino mass- 
eigenstate, if the normal (inverted) mass ordering is re- 
alized, i.e. if Am\ > (Ato^ < 0). The main feature of 
I Ue3 1 is its smallness compared to the other mixing angles 
which makes it difficult to be measured experimentally 
[2]. Three independent groups performing global fits of 
the world's neutrino data find 3a limits of |J7e3p < 0.046 
[3], < 0.053 g], and < 0.043 (0.047) [5], where the two 
values for Ref. stem from two slightly different analy- 
ses of solar data. | Ues \ is of fundamental importance as 
it is crucial to rule out a large number of models con- 
structed to explain the peculiar structure of lepton mix- 
ing [6], CP violation in oscillations depends crucially on 
I Ue3 1 , and it also plays an important role in neutrino-less 
double beta decay [7j. Although some rather weak hints 
on a non-zero value exist [3] at ~1.5f7, all existing data 
are well compatible with |J7e3| = 0. 

The smallness of {Uesl is most naturally explained by 
some flavor symmetry [5HlO) predicting |?7e3| to vanish 
at some high scale A. In general, the symmetry is broken 
and potentially generates small but non-zero values of 
I Ue3 1 at the low scale A relevant for the oscillation exper- 
iments. Radiative corrections are one inevitable source 
of this breaking. In this letter we discuss the minimally 



induced values of |?7e3| through renormalization group 
(RG) evolution, in order to investigate whether non-zero 
values of | C/gs | can be expected in general. Given the huge 
amount of experimental activity and the theoretical im- 
portance of this parameter, it is obviously of interest to 
investigate whether non-zero values can be expected on 
general grounds. We therefore especially focus on the 
case of an inverted mass ordering with = {Uesl ~ Oj 
which is a fixed point at 1-loop. It is shown that at 2-loop 
order this is no longer a fixed point, and we argue that 
it leads to the lowest possible value of |J7e3|, if extreme 
fine-tuning of the parameters is neglected. The order of 
magnitude of \Ue3\ at low scale lies between 10~^^ and 
2Q-i4^ This is, of course, beyond the reach of the planned 
neutrino oscillation facilities, but of interest in supernova 
physics. It is furthermore easy to rule out this possibil- 
ity. We also give the corresponding bounds on \Ue3\ for 
the normal mass-ordering, which happens to be several 
orders of magnitude larger. 

The 1-loop RG evolution of ^13 = sin^^ \ Ue3\, up to 
first order, is given by [11] [T^ 



713 



sin 2^12 sin 26*23" 



m3 



(1) 



327r2 Ami (1-1- i?) 

[mi cos((y3i — 5) — (1 + R)m2 cos((/32 ^ ^ Rm3 cos S] , 

where the dot denotes d/dt = ^d/d^ with t = In {n/A), 
jjL being the renormalization scale. One has C = ~f 
(C ^ 1) in the Standard Model (SM) (Minimal Super- 
symmetric SM (MSSM)), R = b.m\lb.m\, and the tau 
Yukawa coupling is = ^2 (m^/GeV)/246 ~ 0.010 in 
the SM and 0.010 (1 -f tan^ /3) in the MSSM. The elec- 
tron and muon Yukawa couplings ?/e_^ have been ne- 
glected here. Eq. ([T]) shows that the evolution of |?7e3| 
is proportional to ms. Moreover at 1-loop, one has 
rrii oc mi(i = 1,2,3) and thus |[/e3|(Ai) — ?Ti3(/i) — 
is a fixed point of RG evolution: these quantities remain 
zero throughout the RG evolution. The stability of this 
scenario under radiative corrections remains true even if 
?/e.p are taken into account [5]. 
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To understand the existence of the fixed point in a 
more general form, we consider the neutrino mass matrix 
nil, for ma = \Ue3\ = 0, which obeys a structure denoted 
as "scahng" : 



scaline 




(2) 



where tan^ ^23 = ^-nd a, b, d are determined by 

™i,2 5 ^12 J and the Majorana phase difference Aip = ip2 — 
tfi. At 1-loop, the RG equation of is given as pTlfT^ 



1 



167r2 



+ C {Y^ + m,Y)] , (3) 



where Y = Ye,Y} with Yi = diag(?/e, 2//^, J/r) being the 
charged lepton Yukawa matrix. The parameter is 
a function of gauge and Yukawa couphngs (as well as 
the Higgs self-couplings in the SM), which does not lead 
to modifications of the mixing matrix elements. From 
Eq. ([3]) RG evolution of an element (mi,)^^ becomes [T3] 



1 



[a. + C{yl+yl)]im,)^p, (4) 



which at low scale leads to the RG-modified mass matrix 



P 



{m1)er le It 



(5) 



where denotes the mass-matrix at the high-scale and 



/a = exp j^vldt 



(6) 



with a,/3 € {e,/i, r}. Eq. ([s]) shows that the scaling 
property of rrfj'^^™^ is preserved by the RG evolution 
since the second and third column remain proportional 
to each other. The atmospheric neutrino mixing changes 
to tan^ 6*23 = |c^|)- Thus mg = |C/e3| = remains 

valid at all energies during 1-loop running. Therefore, 
any non-zero | Uez \ generated at 2-loop can be considered 
to be the smallest guaranteed value for |C/e3|, which ob- 
viously is of interest in its own right, and should in the 
end represent the final precision goal for any search. If 
this 2-loop value cannot vanish, it follows that in general 
|J7e3| can be expected to be non-zero. 

The largest contribution in Eq. Q will come 
from |C|?/2/(i67r2) ~ 9.5 x 10"^ for the SM and 
~ 6.3 X 10^^ (1-1- tan^ /3) for the MSSM. in Eq. ^ can 
thus be approximated as ~ 1 -f j^-^ In 4 = 1 + . 



At 2-loop, the relevant term in the RG evolution equa- 
tion of nil, in case of the SM is ^3] 



nil 



(167r2 



Y^ m^Y , 



(7) 



H 



H 



T S • 



H 



FIG. 1: Rank-changing 2-loop diagram in the SM which 
breaks the 1-loop RG invariance ms = If^esl = 0. Il {sr) 
are left-(right-)handed leptons and H is the Higgs doublet. 



which arises from the diagram shown in Fig. [TJ In the 
MSSM, there are no such terms (because of the non- 
renormalization theorem) unless supersymmetry (SUSY) 
is broken [15] . Hence in the remaining part of this paper 
we will consider the SM, unless otherwise specified. From 
Eq. ([7| the evolution of (m,j,)^^ can be written as 



I0.P (167r2)2 

and the mass-matrix nin gets modified to 



(8) 



{ml)rr Pr Irr 

(9) 

The new parameters lap are given by 



lap = exp I (ygf^y^ / yl yp dt} 



1 + 



(T6¥5F y/3 in ^ = 1 + eap 



(10) 



As can be seen, 2-loop running increases the rank of 
from 2 to 3. In Eq. ( [To| the largest contribution comes 
from I err I — 2.0 x 10^'"^, and thus numerically is of the 
same order as that of at 1-loop. Here we have taken 
A = 102 GcV and A = 10^2 GeV, which wiU be used 
throughout this paper. Obviously, the scale of e,-r sets 
the value for |J7e3| generated at 2-loop order. 

To explicitly calculate the implied smallest value of 
jC/eal after 2-loop running, we consider the scaling sce- 
nario 7713 = iJ/esI = at the high scale A. It suffices 
to take only e^-^. into account and neglect all other tap 
in Eq. (10). As a result of the change in rank of m^. 



the 2-loop running, unlike the 1-loop case, gives rise to a 
non-zero smallest neutrino mass 



ma 



\ siv? 2023 {ni2 cos2 6*12 e 



i m2 sin^ 2023 (cos2 612 e 



(1- 



- sin 6*12 mi) {errl 
f) sin2 0,2) \err\ 

-13 



The order of ma is ~ m2 |erT|/4 ^ 2.5 x lO"^"* eV for 
m2 ~ 0.05 eV. It confirms the result from Ref. [T3]. 

The 2-loop evolution equation for |C/e3| in the scaling 
scenario can be obtained from Eq. ([7| , following the pro- 



3 




FIG. 2: Generated value of \Ues\ at A = 100 GeV assuming a 
vanishing value at A = 10^^ GeV in the leading-log approxi- 
mation. The (solid) red line corresponds to fixing all mixing 
parameters to the best-fit values and the yellow region shows 
the 3(7 region, which are taken from [4]. 



cedure of [16] (see also [l7|), to be 



dt 



e3| 



\Ue: 



Vt 



'23 



2 (IBtt^)^ mi 7712 



sin 26*12 sin 26*23 X 



mi — m,2 e 



I "^1 ^12 + 1^2 Ci2 e.^^^ I 



in the limit ye,^i Vt- The /3-function of \Ue3 
depend on the Dirac CP phase S for vanishing 
which has been discussed in detail in [T5] 
the typical order of magnitude of \Ue3 



(11) 
does not 

Ue3\ m, 

From Eq. ( [TT| ), 
at low scale A 



is \Ue3\ ^ 2/^/(4 (167r2)2) In A/A ~ IQ-^^. Fig. § shows 
the dependence of \Ue3\ on the Majorana phase difference 
A(p. The maximal value of |J7e3| is obtained for cos Aip — 
(tan^ 012 — l)(m| — mi tan^ 6'i2)/(4mi m2 tan^ ^12), and 
gives 



\U, 



e3\ 



3.2 X 10" 



(12) 



while the minimum value, obtained for A(p = 0, is 



\U, 



e3\ 



2.0 X 10" 



It is suppressed with respect to Eq. (12 1 by a factor of 
order R ~ 1/30, which remains true even if ?/e,p are taken 
into account. It is clear that a 3-loop effect could not 
cancel this suppressed value, since it would be further 
suppressed by l/(167r^) <^ R. 

One possibility for a smaller value of \Ue3\ would be if 
an initially non-zero value accidentally runs to (almost) 
zero in the course of its evolution which, of course, will 
need an extreme amount of fine-tuning. Another possibil- 
ity arises when the high energy completion of the effective 
theory is considered and the RG effects in the complete 
theory conspire with the effective running to give tiny 



values, e.g. threshold effects in case of seesaw models. 
However, that will depend on the particular high energy 
extension considered, and will again involve extreme fine- 
tuning. For example, in seesaw models threshold effects 
generally lead to larger RG running than the running in 
the effective theory [T7] , and hence a fine-tuned cancella- 
tion will be required to produce the tiny \Ue3 \ of relevance 
in this discussion. Hence it is natural to claim that the 



2-loop generated |i7e3l as given in Eq. ( 11 ) is the smallest 



value one should expect. Moreover, it shows that non- 
zero \Ue3\ can be expected on general grounds. 

As previously mentioned, \Ue3\ is not generated at 2- 
loop in the MSSM, as long as SUSY is unbroken. How- 
ever, analogous to the discussion for ma in [l4], |J7e3| is 
generated as SUSY is broken and receives contributions 
of two different origins. First, there is a 2-loop contri- 
bution with sleptons and selectrons in the loop as shown 
in Fig. 1 of [14 . Although it is not enhanced by a large 
logarithm, it depends on (1 -|-tan^ /3)^ along with a com- 
plicated order one function of sparticle masses. Second, 
there can be corrections due to the the SUSY breaking 
operator LHuLHu [19 , which leads to \Ue3\ of similar 
order as the 2-loop contribution. Thus while a general 
prediction is not possible, it suffices to estimate the order 
of magnitude to be 



It/, 



e3|™(l + tan2/3)Vlny, 



(14) 



SM corresponds to the values in the SM dis- 

MSSM 



where \Ue3 

cussed above. Thus even for moderate tan/3, \Ue3, 
will be larger than the corresponding SM value. 

Can such values be tested? With the present view, the 
ultimate experiment in order to study neutrino oscilla- 
tions will be a neutrino factory. However, the facilities 
currently under study [T] have a 3(7 discovery potential 
on \Ue3\ of at most 1.5 x 10"'^, and hence the 2-loop gen- 
erated tiny \Ue3\ is far beyond its scope. 

However, tiny non-zero values of |C/e3|j as discussed 
here, may have observable consequences in the neutrino 
spectra emitted by supernovae. Extremely high neutrino 
densities around the neutrinosphere give rise to so-called 
collective effects ^ , of which "spectral swapping" is the 
(13) one where Vg and v^ ^ swap their energy spectra at high 
(low) energies but keep their spectra at low (high) ener- 
gies in case of inverted (normal) hierarchy. Here "high" 
and "low" energies are understood relative to a critical 
energy Ec, which decreases with \Ue3\ in the normal or- 
dering, while is independent of \Ue3\ if the mass-ordering 
is inverted [23]. The effect can be well understood with 
the analogy of an inverted pendulum [25] . which turns 
around by the slightest instability, generated by a non- 
zero \Ue3\ in this case. In case of inverted mass-ordering, 
the spectral swap occurs regardless of how small but non- 
zero |C/e3| is, the impact of |t/e3| being only to introduce 
a logarithmic dependence on the radius at which the con- 
version sets in. Ref. |MI discusses in this context values 
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of lUesl down to 10"^", 56 orders of magnitude smaller 
than the lower limit obtained here. 

While the discussion here may be a bit too simplified 
(the effects may depend on supernova details, they can 
be also induced by other effects like an asymmetry in the 
primary v^-Vt fluxes as well as radiatively induced mat- 
ter effects for v^^r etc.), it catches the main point to 
be made here: tiny values of |J7e3|, in particular the ones 
generated by 2-loop running in the inverted hierarchy, are 
expected to have observable consequences in supernova 
physics. The lower bounds on |J7e3| presented here may 
be helpful in analyses of supernova neutrino analyses. 

The scenario leading to the smallest possible value of 
|C/e3| can easily be tested or ruled out, not only by long- 
baseline oscillation experiments pinning down the mass 
ordering, but also via neutrino mass related observables: 
for ma = |J7e3| = the effective mass governing neutrino- 
less double beta decay is simply 



Ami I 



1-sin^ 26*12 sin2Av3/2. 



(15) 



The Majorana phase difference A(/3 is the same as the 
one determining the magnitude of |C43| and is, in princi- 
ple, measurable by precise measurements of double beta 
decay [20]. In contrast to its value in the normal or- 
dering, (m) is bounded from below by -y/jAmJI cos26'i2, 
which is testable in future experiments |21j . In what re- 
gards cosmology, the sum of neutrino masses is about 
2-\/|Ami| ~ 0.1 eV, which is twice as high as its value 
in the normal ordering, and might be testable in future 
measurements as well [22] . 

In case of normal ordering the minimal | C/e3 1 can triv- 
ially be obtained from the 1-loop running given in Eq. ^. 
The minimal value is obtained with Lp i — Lp2 = to be 



\Ue3\ = Sin 2012 sin 2023- 



m3 (m2 - mi) A 
' (mi -f ma) (m2 -f ma) A 
(16) 

which quantifies to \Ue'i\ ^ 1.1 x 10 ^ in the SM and 
~ 1.6 X 10"^ (1 + tan2/3) in the MSSM. Interestingly, 
jC/eal is smallest for quasi-degenerate neutrinos and not 
for mi = 0. Introducing the common mass scale mo 
of the quasi-degenerate neutrinos, the lowest values are 
decreased by a factor of 



(/ 



AmQ 



1600 



AmQ 



vo.sevy 



and hence reduce the scale of \Uf,^\ by at most three 
orders of magnitude, and still stay three orders above 
the 2-loop value in the inverted hierarchy. 

In summary, we have analyzed the 2-loop running of 
Ue3 in the effective low energy theory and showed that it 
gives the minimal value of |J7e3| in scenarios of inverted 
mass hierarchy with Uez = at high scale. We focussed 



on this scenario because at 1-loop there is no generation 
of a non-zero value, and because it is of interest to 
investigate whether non-zero \Ued,\ can be expected in 
general. The order of magnitude of jC/esI generated in 
this way is between 10~^^ and lO"^'* depending on the 
Majorana phases. Though such small values are beyond 
the reach of currently foreseen neutrino oscillation 
experiments, they are obviously of fundamental interest, 
and possess applications in supernova physics. 
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